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Intisari 

Proses difusi menjadi salah satu mekanisme penting di dalam sistem tubuh manusia, seperti pada sistem saraf 
yaitu difusi Ca 2+ yang terjadi pada sinapsis. Tujuan penelitian ini adalah mempelajari proses difusi Ca 2+ 
yang terjadi jika Gd menghalangi saluran Ca 2+ pada pra-sinapsis. Dalam penelitian ini tidak menggunakan 
Gd tetapi dimodelkan dengan jumlah saluran Ca 2+ yang bervariasi, semakin sedikitnya jumlah saluran Ca 2+ 
yang terbuka menggambarkan semakin banyak Gd yang menghalangi saluran Ca 2+ . Simulasi ini menggu¬ 
nakan metode Monte Carlo Cell (MCell). Model sinapsis dibuat dengan Blender 2.76b yang terdiri dari dua 
kompartemen yang digambarkan dengan dua lingkaran. Molekul di kompartemen luar pra-sinapsis berinter- 
aksi dengan molekul permukaan kompartemen pra-sinapsis dan menghasilkan molekul di dalam pra-sinapsis. 
Digunakan tiga indikator dalam simulasi untuk mengetahui pengaruh proses difusi yang terjadi pada sinapsis, 
hasil simulasi untuk molekul mencapai jumlah setengah pada variasi jumlah luasan saluran dari kecil ke besar 
berturut-turut membutuhkan waktu 0,65 s, 0,34 s, 0,24 s, dan 0,18 s. Waktu untuk variasi kerapatan dari besar 
ke kecil berturut-turut 0,337 s, 0,362 s, dan 0,364 s. Jumlah molekul terdifusi pada variasi backward rate dari 
besar ke kecil dengan waktu yang sama yaitu 3928, 3626, 2037, dan 529 molekul. 

Abstract 

The diffusion process became one of the important mechanisms in the human body system, such as in the 
nervous system which is the Ca 2+ diffusion in the synapses. This research aims to study the Ca 2+ diffusion 
process if Gd inhibits Ca 2+ channel in pre-synapses. This study did not use Gd, but it was modeled with 
various numbers of Ca 2+ channels. The less number of Ca 2+ opened channels described the more Gd which 
inhibited Ca 2+ channels. The Monte Carlo Cell Method (MCell) was used in this study. The synapses model 
was made using Blender 2.76b; consisted of two compartments illustrated with two circles. The molecule in the 
outside pre-synapses compartment interacted with the molecular surface of the pre-synapses compartment and 
resulted in a molecule in pre-synapses. The result of the simulation process for the molecules reached half of 
the variations in the number of total area channels. The time needed being 0.65s, 0.34s, 0.24s, and 0.18s, from 
the smallest to the biggest respectively. The times needed for density variations from big to small were 0.337s, 
0.362s, and 0.364s respectively. The number of the diffused molecules in the backward rate variation within the 
same time was 3928, 3626, 2037, and 529 molecules respectively. 
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I. PENDAHULUAN 


Penggunaan Magnetic Resonance Imaging (MRI) di dunia 
kesehatan semakin berkembang, secara umum penggunaan 
MRI ini membutuhkan kontras agen yaitu Gadolinium (Gd). 
Gd masuk ke dalam tubuh dengan cara disuntikkan ke dalam 
tubuh dan terdifusi ke bagian-bagian tubuh. Gd memiliki ke- 
untungan dapat meningkatkan kualitas gambar MRI yang di- 
hasilkan [1]. Meskipun demikian, dalam 10 tahun terakhir 
ini ditemukan bahwa Gd yang diberikan pada tubuh tidak 
sepenuhnya keluar melalui proses eksresi tubuh tetapi ada 
Gd yang tertinggal di dalam otak. Otak mengandung banyak 


neuron yang memiliki kemampuan masing-masing, salah satu 
proses difusi Gd yang tertinggal dalam otak terjadi pada sinap¬ 
sis [2]. 

Sinapsis adalah titik temu antara terminal akson salah satu 
neuron dengan terminal dendrit neuron lain yang memung- 
kinkan pertukaran informasi antar neuron. Terdapat dua 
sinapsis di dalam tubuh manusia yaitu kimia dan elektrik. 
Pembahasan ini hanya fokus pada sinapsis kimia karena difusi 
banyak terjadi pada sinapsis kimia. Sinapsis kimia melakukan 
komunikasi melalui sekresi neurotransmiter dengan mengak- 
tifkan sel-sel. Sistem sinapsis terdapat dua bagian penting 
yaitu pra-sinapsis dan pascasinapsis. Pada permukaan pra- 
sinapsis terdapat saluran yang disebut sebagai saluran Ca 2+ . 
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Gambar 1: (a). Ilustrasi sinapsis dan saluran Ca 2+ sebagai tempat terjadinya difusi [5], (b). Model pra-sinapsis yang dibuat 

menggunakan Blender. 


Saluran Ca 2+ ini berfungsi sebagai pintu masuknya Ca 2+ dan 
akan terbuka jika adanya potensial aksi [3]. Ca 2+ masuk dari 
luar pra-sinapsis ke dalam pra-sinapsis secara difusi sehingga 
menimbulkan peningkatan kadar Ca 2+ di dalam terminal pra- 
sinapsis [4]. Jumlah saluran Ca 2+ yang terbuka setiap waktu 
memiliki korelasi dengan jumlah saluran Ca 2+ yang tidak ter- 
halang oleh suatu materi lain (dalam kasus ini adalah Gd). Gd 
menghalangi saluran Ca 2+ maka jumlah saluran Ca 2+ yang 
terbuka semakin berkurang [5,6]. 

Proses difusi Ca 2+ yang terjadi pada sinapsis ini menjadi 
salah satu mekanisme penting di dalam sistem tubuh manusia. 
Difusi ini terjadi karena adanya perbedaan jumlah molekul an- 
tara di luar pra-sinapsis dan di dalam pra-sinapsis sehingga 
molekul Ca 2+ berpindah dari konsentrasi tinggi ke konsen- 
trasi rendah. Proses difusi berakhir ketika konsentrasi di luar 
pra-sinapsis dan di dalam pra-sinapsis berada pada keadaan 
setimbang [4,7]. Peristiwa difusi yang terjadi dapat diny- 
atakan dengan Hukum Fick yaitu, pada arah tertentu molekul 
Ca 2+ yang melewati suatu luasan tertentu tiap unit satuan 
waktu adalah sebanding dengan gradien konsentrasi pada arah 
tersebut [8]. Molekul Ca 2+ yang melewati luas saluran Ca 2+ 
tiap unit satuan waktu dapat dinyatakan dalam Pers.(l) [9]: 

bC 

( 1 ) 

Penelitian ini bertujuan untuk mempelajari proses difusi 
Ca 2+ yang terjadi jika Gd menghalangi saluran Ca 2+ pada 
pra-sinapsis, dalam penelitian ini tidak menggunakan Gd 
tetapi dimodelkan dengan jumlah saluran Ca 2+ yang bervari- 
asi, semakin sedikitnya jumlah saluran Ca 2+ yang terbuka 


menggambarkan semakin banyak Gd yang menghalangi salu¬ 
ran Ca 2+ . Simulasi yang digunakan adalah metode Monte 
Carlo Cell (MCell), yang dapat mensimulasikan gerakan dan 
reaksi molekul di dalam dan di antara sel-sel berdasarkan 
proses difusi yang dikombinasikan dengan proses biolo- 
gis, MCell dapat juga menghitung semua probabilitas yang 
mungkin terjadi dengan hasil yang presisi dikarenakan adanya 
batas secara geometri [10]. 


II. METODOLOGI 

Simulasi MCell dilakukan menggunakan MCell 3.3 yang 
dijalankan dengan dengan PC (.Personal Computer) Win¬ 
dows 7 Intel Core i7-6700 CPU @ 3,40 GHz GPU NVIDIA 
GeForce GTX 960. Visualisasi sinapsis (luar pra-sinapsis 
dan dalam pra-sinapsis) dibuat dengan menggunakan Blender 
2.76b yang add-on dengan Cellblender 1.0. 

Model pra-sinapsis yang dapat diilustrasikan pada Gam- 
bar 1(a). Dalam simulasi ini didekati dengan model ico 
sphere yang dibuat dengan Blender 2.76b yang terdiri dari 
dua kompartemen yang digambarkan dengan dua lingkaran 
yang terdapat pada Gambar 1(b). Lingkaran besar adalah 
luar pra-sinapsis dengan volum 1,02 /im 3 dan luas permukaan 
4,90 ptm 2 . Lingkaran kecil adalah dalam pra-sinapsis dengan 
volum 0,11 /im 3 dan luas permukaan 1,13 ptm 2 . Karakter- 
istik dasar dinding permukaaan pra-sinapsis adalah reflektor 
atau pemantul sehingga dibuat dinding model permeabel de¬ 
ngan suatu luasan. Lubang pada permukaan berfungsi sebagai 
saluran Ca 2+ sehingga molekul Ca 2+ di luar pra-sinapsis bisa 
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Gambar 2: Hasil simulasi menggunakan MCell. Ca 2+ masuk 
dari luar pra-sinapsis (garis hitam) ke dalam pra-sinapsis 
(garis merah) melalui saluran Ca 2+ . 


masuk ke dalam pra-sinapsis. 

Molekul di kompartemen luar pra-sinapsis (A) berinter- 
aksi dengan molekul permukaan kompartemen pra-sinapsis 
(ti) dan menghasilkan molekul di kompartemen dalam pra- 
sinapsis (B). Simulasi ini menggunakan proses reversible 
yang memungkinkan terjadinya perpindahan molekul Ca 2+ 
dari dalam pra-sinapsis ke luar pra-sinapsis. 

A-\-t\ B (2) 

Digunakan tiga indikator dalam simulasi untuk mengetahui 
proses difusi yang terjadi pada sinapsis. Simulasi pertama 
adalah dengan memvariasikan jumlah saluran Ca 2+ yaitu 8, 
16, 24 dan 32, saluran ini berbentuk segi enam dengan re- 
septor terdiri dari 6 element grip (EG) yang berbentuk segit- 
iga [5], sehingga jumlah salurannya menjadi 48, 96, 144, dan 
192 segitiga dengan luasan tiap segitiga adalah 2,2 x 10 -4 
/im 2 . Molekul permukaan saluran Ca 2+ yang digunakan yaitu 
10, 20, 30, dan 40 agar memiliki nilai kerapatan yang sama, 
dan molekul representatif dari Ca 2+ yang berada di luar pra- 
sinapsis (molekul volum) menggunakan jumlah yang tetap 
yaitu sebesar 5000 molekul. 

Simulasi kedua adalah dengan memvariasikan besar kerap¬ 
atan saluran Ca 2+ dengan mengubah-ubah jumlah saluran 48, 
96, dan 192 dan menggunakan molekul permukaan dengan 
jumlah yang tetap yaitu sebesar 10 molekul sehingga jika di- 
hitung besar kerapatannya akan memiliki nilai yang berbeda. 
Molekul volum digunakan jumlah yang tetap yaitu sebesar 
5000 molekul. 

Simulasi ketiga adalah dengan memvariasikan nilai back¬ 
ward rate yaitu 10 5 , 10 6 , 10 7 , dan 10 8 M -1 s -1 . Jumlah salu¬ 
ran Ca 2+ yang digunakan adalah 144 saluran, molekul per¬ 
mukaan sebesar 30 molekul dan molekul volum sebesar 5000 
molekul. 


III. HASIL DAN PEMBAHASAN 

Gambar 2 menunjukkan distribusi molekul Ca 2+ sebe¬ 
sar 5000 molekul dari luar pra-sinapsis (garis hitam) ke 
dalam pra-sinapsis (garis merah) melalui saluran Ca 2+ de¬ 
ngan jumlah saluran sebanyak 144 dan molekul permukaan 
sebesar 30 molekul. Jumlah molekul Ca 2+ di luar pra- 
sinapsis menurun secara bertahap karena molekul berdifusi 
ke dalam pra-sinapsis sehingga terjadi penambahan molekul 
Ca 2+ di dalam pra-sinapsis secara bertahap. Pada saat t 
= 0.24 s, molekul Ca 2+ di luar dan di dalam pra-sinapsis 
memiliki nilai yang sama yaitu sebesar 2500 molekul. Sete- 
lah sampai pada keadaan setimbang tersebut pertukaran Ca 2+ 
dari luar ke dalam pra-sinapsis atau sebaliknya masih terjadi 
karena adanya dorongan dengan menggunakan forward rate 
dan backward rate yang berbeda yaitu 10 8 M -1 s -1 dan 10 6 
M -1 s -1 . 

A. Sutresno, et al. [4], menyatakan bahwa metode Monte 
Carlo berhasil mensimulasikan interaksi antara Ca 2+ dan Gd 
sebagai proses difusi di sinapsis. Hasil penelitian tersebut 
menyediakan kerangka kerja untuk menentukan rasio tertentu 
antara Ca 2+ dan Gd sehingga proses difusi masih dapat ter¬ 
jadi. Kecepatan difusi Ca 2+ ke presinaps dan kecepatan 
Gd memblokir saluran Ca 2+ tergantung pada rasio keduanya. 
Penelitian ini merepresentasikan hal tersebut, yaitu pengaruh 
saluran Ca 2+ yang terblokir oleh Gd dalam proses difusi 
di pra-sinapsis tetapi dalam hal ini tidak menggunakan Gd. 
Telah didapat hasil simulasi untuk pengaruh jumlah saluran 
Ca 2+ , pengaruh kerapatan saluran Ca 2+ , serta pengaruh for¬ 
ward rate dan backward rate dalam proses difusi yang terjadi 
di pra-sinapsis sebagai berikut: 

Pengaruh Jumlah Saluran Ca 2+ 

Pengaruh jumlah saluran Ca 2+ terhadap proses difusi dike- 
tahui dengan menggunakan simulasi variasi jumlah lubang 
saluran Ca 2+ 48, 96, 144 dan 192 saluran, kerapatan dibuat 
tetap dengan menggunakan molekul permukaan berturut-turut 
10, 20, 30, dan 40 molekul. Variasi jumlah saluran Ca 2+ 
menginterpretasikan Gd yang menghalangi saluran Ca 2+ . 
Gambar 3 menunjukkan distribusi 5000 molekul Ca 2+ dengan 
waktu yang dibutuhkan molekul untuk berdifusi. 

Grafik yang terdapat pada Gambar 3 menunjukkan bahwa 
dengan jumlah molekul yang sama dan jumlah saluran 
yang berbeda dapat mempengaruhi waktu yang dibutuhkan 
molekul untuk berpindah ke dalam pra-sinapsis. Jumlah salu¬ 
ran yang berbeda berarti luasan saluran juga berbeda karena 
luas saluran diperoleh dengan mengalikan jumlah saluran de¬ 
ngan luasan setiap saluran. Pada tiap luasan saluran yang 
berbeda, molekul mencapai jumlah dari setengahnya terjadi 
pada waktu yang berbeda. Untuk mencapai 2500 molekul, 
jumlah saluran 48 membutuhkan waktu 0,65 s, saluran 96 
membutuhkan waktu 0,34 s, saluran 144 membutuhkan waktu 
0,24 s, dan saluran 192 membutuhkan waktu 0,18 s. Dari hasil 
simulasi tersebut menunjukkan berkurang dan bertambahnya 
jumlah saluran berpengaruh pada kecepatan difusi Ca 2+ , se¬ 
hingga semakin sedikit jumlah saluran Ca 2+ yang terbuka 
maka semakin lambat Ca 2+ berpindah ke dalam pra-sinapsis. 
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Gambar 3: Jumlah molekul terhadap waktu dengan lubang 
48 (garis hitam), lubang 96 (garis merah), lubang 144 (garis 
biru), dan lubang 192 (garis hijau). 



Vs < detik ) 

Gambar 4: Kerapatan saluran Ca 2+ terhadap to, 5 . 


Jumlah saluran yang berkurang menggambarkan jumlah salu¬ 
ran di pra-sinapsis yang tertutup oleh Gd dan membuat difusi 
Ca 2+ terganggu karena luas saluran Ca 2+ lebih kecil. 

Pengaruh Kerapatan Saluran Ca 2+ 

Pengaruh kerapatan saluran Ca 2+ terhadap proses difusi dike- 
tahui dengan menggunakan jumlah saluran 48, 96, dan 192, 
molekul volum sebesar 5000 molekul, dan menggunakan 
molekul permukaan dengan jumlah yang tetap yaitu sebesar 
10 molekul. Besarnya kerapatan diperoleh dengan menghi- 
tung jumlah molekul dibagi luas permukaan sehingga be- 
sar kerapatan berturut-turut adalah 5700, 2800, dan 1400 
fim 2 . Gambar 4 menunjukkan hubungan kerapatan saluran 
dengan to, 5 . Grafik yang terdapat pada Gambar 4 menun¬ 
jukkan bahwa dengan nilai kerapatan yang berbeda mem- 
pengaruhi waktu yang dibutuhkan molekul volum mencapai 



Gambar 5: Jumlah molekul terhadap waktu untuk backward 
rate 10 5 M - 1 s -1 (garis hitam), 10 6 M - 1 s -1 (garis merah), 
10 7 M - 1 s -1 (garis biru), dan 10 8 M - 1 s -1 (garis hijau). 


jumlah setengahnya. Kerapatan berhubungan dengan jum¬ 
lah molekul yang ada di saluran Ca 2+ , semakin besar kera- 
patannya maka peluang molekul volum untuk berinteraksi de¬ 
ngan molekul permukaan semakin besar sehingga waktu yang 
dibutuhkan molekul untuk berpindah lebih cepat dan molekul 
volum yang masuk lebih banyak. Untuk mencapai 2500 
molekul, kerapatan 5700 fim 2 membutuhkan waktu 0,337 s, 
kerapatan 2800 /im 2 membutuhkan waktu 0,362 s, dan ker¬ 
apatan 1400 /im 2 membutuhkan waktu 0,364 s. Hasil simu- 
lasi tersebut menunjukkan semakin kecil kerapatannya, waktu 
yang dibutuhkan lebih lama. Kerapatan yang semakin berku¬ 
rang menggambarkan molekul permukaan di pra-sinapsis ter- 
halang oleh Gd dan membuat difusi Ca 2+ terganggu karena 
probabilitas interaksi antara molekul volum dan molekul per¬ 
mukaan lebih kecil. 

Pengaruh Forward Rate dan Backward Rate 

Pengaruh backward rate terhadap proses difusi diketahui de¬ 
ngan menggunakan simulasi variasi nilai backward rate 10 5 , 
10 6 , 10 7 , dan 10 8 M - 1 s -1 , dengan jumlah saluran sebanyak 
144, forward rate 10 8 M - 1 s -1 , dan molekul volum sebesar 
5000. Gambar 5 menunjukkan distribusi 5000 molekul Ca 2+ 
dengan waktu yang dibutuhkan molekul untuk berdifusi. 

Forward rate dan backward rate menunjukkan ke- 
cepatan reaksi dari luar pra-sinapsis menuju dalam pra- 
sinapsis, semakin besar perbandingan perbedaannya maka 
proses difusinya lebih cepat. Kemiringan dari setiap grafik 
menunjukkan perbedaan jumlah molekul dengan waktu yang 
berhubungan dengan fluks yaitu molekul Ca 2+ yang melewati 
luas saluran Ca 2+ tiap unit satuan. Pada waktu 0,5 s untuk 
setiap variasi backward rate menghasilkan jumlah molekul 
yang berbeda, pada backward rate 10 5 M - 1 s -1 molekul yang 
terdifusi mencapai 3928 molekul, pada backward rate 10 6 
M - 1 s -1 molekul yang terdifusi mencapai 3626 molekul, pada 
backward rate 107 M - 1 s -1 molekul yang terdifusi mencapai 
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2037 dan pada backward rate 10 8 M -1 s -1 molekul yang ter- 
difusi hanya mencapai 529 molekul. Hasil simulasi tersebut 
menunjukkan dengan besar waktu yang sama semakin besar 
backward rate, fluksnya semakin kecil karena molekul volum 
yang berdifusi semakin sedikit. 

IV. SIMPULAN 

Metode MCell dapat digunakan untuk mensimulasikan 
proses difusi Ca 2+ pada sinapsis sehingga dapat mempelajari 
proses difusi yang terjadi. Jumlah saluran Ca 2+ , kerapatan 
saluran Ca 2+ dan backward rate mempengaruhi proses difusi 
Ca 2+ pada pra-sinapsis. Jumlah saluran Ca 2+ yang berbeda 
48, 96, 144, dan 192 untuk mencapai jumlah molekul yang 
sama memiliki waktu yang berbeda secara berturut-turut yaitu 
0,65 s, 0,34 s, 0,24 s, dan 0,18 s. Semakin banyak jumlah salu¬ 


ran Ca 2+ berarti semakin luas saluran Ca 2+ maka waktu yang 
dibutuhkan molekul untuk berpindah semakin cep at. Besar 
kerapatan saluran Ca 2+ dengan jumlah molekul volum tetap 
berbanding lurus dengan laju difusi, untuk mencapai 2500 
molekul dengan kerapatan 5700 ptm 1 2 membutuhkan waktu 
0,337 s, kerapatan 2800 fi m 2 4 membutuhkan waktu 0,362 s, 
dan kerapatan 1400 ptm 2 membutuhkan waktu 0,364 s, se¬ 
hingga semakin besar kerapatan saluran maka semakin cepat 
perpindahan molekul Ca 2+ . Perbandingan perbedaan for¬ 
ward rate dan backward rate mempengaruhi waktu yang dibu¬ 
tuhkan molekul untuk berpindah,saat waktunya 0,5 s pada 
backward rate 10 5 M -1 s -1 mencapai 3928 molekul, pada 
backward rate 10 6 M -1 s -1 mencapai 3626 molekul, pada 
backward rate 10 7 M -1 s -1 mencapai 2037 dan pada back¬ 
ward rate 10 8 M -1 s -1 hanya mencapai 529 molekul, se¬ 
hingga semakin besar perbandingannya maka molekul yang 
berpindah semakin sedikit. 
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